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Abstract: We investigate the properties of holographic fermions in charged Lifshitz 
black holes at finite temperature through the AdS/CFT correspondence. In the charged 
Lifshitz background with the dynamical exponent z = 2, we find that the dispersion 
relation is linear. The scaling behavior of the imaginal part of the Green function 
relative to = k — k F is also discussed. We find, although the system has linear 
dispersion relation and quadratic quasi-particle width, it does not satisfy Luttinger's 
theorem. We also find that the variation of the scaling parameters a and f3 is small as 
the charge q varies. Furthermore, we also discuss the effect of the dynamical exponent 
z by considering the cases z = 4 and z = 6 and show that ImGu become smooth when 
the dynamical exponent z increases. 
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1. Introduction 

Recently, there has been an increased interests in studying strongly interacting fermionic 
systems by a gravity dual using the holographic Anti-de Sitter /conformal field theory 
correspondence (AdS/CFT) or Anti-de Sitter/condensed matter theory (AdS/CMT) 
correspondence ^J, |], [3[]. In [|J [5], ||, [7|, it was shown that holographic duality can de- 
scribe Fermi and non- Fermi liquids. In general, the non-Fermi liquid has a sharp Fermi 
surface but its gapless charged excitations differ significantly from those predicted from 
Landau's Fermi liquid theory. Up to now, a proper theoretical framework characteriz- 
ing non- Fermi liquids is lacking. Fortunately, the AdS/CMT correspondence provides 
a large theoretical playground of non-Fermi liquids and give a handle on a difficult 
strong coupling problem. 

In many condensed matter systems, one can find phase transitions governed by fixed 
points with Lifshitz dynamical scaling t — > X z t, x — > Xx. Most interesting examples 
arising in the condensed matter systems correspond to the anisotropic scale invariance 
with z = 2 that describes multicritical points in certain magnetic materials and liquid 
crystals, although when z = 1 the scaling describes scale invariance which arises in the 
AdS/CFT correspondence. In fact, there are many non-relativistic z ^ 1 fixed points 
in condensed matter system including ulta-cold atom mixtures, Fermions at unitarity 
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and Bi-layer graphene. The gravity dual of such non-relativistic field theories was first 
realized by Kachru, Liu, and Mulligan in ||. 

In this paper, we study the fermionic systems numerically in charged Lifshitz theory 
at finite temperatures. Black holes in asymptotically Lifshitz spacetime provide a 
window onto finite temperature effects in strongly coupled Lifshitz models. Black hole 
solutions and their related properties in Lifshitz spacetime were studied by several 
authors || 10, 12, 13, 14, 15, 16, 17 . Charged Lifshitz black hole solutions were 
constructed in fl8| by adding a second Maxwell field J-m) and a scalar field ip which 
was charged under the new gauge field but neutral under the original gauge fields F( 2 ) 
and H(2) ■ Later, this was generalized to an arbitrary (d + 2) dimensional spacetime 



in [ 19[ . A much more general class of solutions of charged Lifshitz black holes was 



discussed numerically in |20| . In the very recent paper , the retarded single-particle 
Green's function is derived in the Lifshitz background for a fermion coupled to the 
conformal field theory at vanishing spatial momenta, temperature and in the absence 
of the chemical potential. The characters of non-relativistic fermionic retarded Green's 
function in Lifshitz geometry with critical exponent z = 2 [|Tj| were studied in |Z| . In 
this paper, we mainly focus on the charged Lifshitz black hole and its gravity duals 

2 (i.e. three-dimensional black holes), z = 4 
i.e. five- dimensional black holes). We 



19] with dynamical exponents z 



given m 

(i.e. four-dimensional black holes) and z 



6 



find that the charge q and the exponents z can affect the properties of the holographic 
fermions in different ways. 

The rest of the paper is organized as follows: We begin by giving a brief review of the 
charged Lifshitz black hole solutions. In section 3, we investigate the Dirac equation as 
well as the boundary conditions in charged Lifshitz black hole backgrounds. By solving 
the Dirac equation numerically, we discuss the properties of the holographic fermions, 
such as the values of the Fermi momenta kp and scaling behaviors with different charge 
q, in z = 2 charged Lifshitz black hole background in the first part of section 4. After 
that, we turn to z = 4 and z = 6 cases. We present the conclusion in the last section. 



2. Charged Lifshitz Black Holes 

In this section, we briefly review on the charged Lifshitz black hole solutions. We start 
with the d + 2 dimensional spacetime with the action 

S = y^t^ / d d+2 xV=9~(R - 2A - \f, v F^ - l -ml ector A,A» - (2-1) 

where F^ u and J-" M „ are gauge fields strength. Here is introduced to obtain charged 
black hole solution under the massive "gauge field" in this spacetime. Meanwhile 
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F^ v plays a key role in modifying the asymptotic geometry from AdS to Lifshitz. Note 
that the massive vector field cannot directly be associated with a global £7(1) gauge 
field because there are no such gauge degrees of freedom, though its kinetic terms look 



like a Maxwell term [IS, |19j. Therefore, we will utilize the J-^ field with global £7(1) 
symmetry in which a non-zero charge density can be introduced in what follows. 
Then based on the action fl2.1| ), the black hole solution has the form[TE|, [T9| 



„2 d 



/(r) = l-4 = l~ Ql 



r z 2d 2 r z 



A t = ^^-^Lr z f(r), F rt — QLr z ~ d ~ x , m 2 vector = ^. (2.2) 

which is asymptotical to the Lifshitz spacetime as r goes to infinity. Here z and Q 
are dynamical exponent and charge parameter, respectively. r + denotes the position 
of horizon which satisfies f(r + ) =0. As discussed in it is worth noting that the 
solution ( |2.2j ) is only valid for z = 2d > 1. We emphasize that the charged Lifshitz 
black hole is very different from the charged black holes in AdS spacetime because it 
does not admit extremal solutions. Thus the behavior of charged Lifshitz black holes 
looks more analogous to Schwarzschild-AdS black holes than Reissner-Nordstrom-AdS 
(RN-AdS) black holes. From in the above solution (|2.2| ), we can obtain the gauge 
field 

A = ^(1 - %). (2.3) 

Note that when z — 1 and d = 2, (2.3) reduces to the scalar potential of charged black 
holes. We stress that the vector field A t is massless and should be used in the following 
discussions of holographic fermions. The metric is invariant under the Lifshitz scaling 

t -)> r z t, I -)> Ir, Xi -)> r 1 ^. (2.4) 



The Hawking temperature at the event horizon is given by 

z 

47T 

where 



T = — r+, (2-5) 



(2.6) 



, _ 
T+ 2d 2 ' 

The Hawking temperature is also invariant under the transformation T — > l z T. The 
entropy of the black hole is 

Sbh = (2.7) 



4G, 



cZ+2 
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with Vd the volume of the d- dimensional spatial part. 

In the rest of this work, we will apply the black hole solution described above to 
discuss some properties of the holographic fermions in this background. Although the 
charged Lifshitz black hole has no extremal limit, we expect that at a low but finite 
temperature we can find evidence of Fermi surfaces. 



3. Dirac Equation 

In order to study fermions in the dual boundary theory, we consider the bulk action 
for a probe Dirac fermion with the mass m and charge q. 1 The bulk fermion action is 



S bu ik = J d d+2 x^i^{T a D a 



m)ip 



(3.1) 



where T a = (e„) a r M and the covariant derivative D a is 



D a = d a + -(u^) a Y^ -iqA ai 



(3.2) 



with = §[1^,1^]. Here, a, b are usual spacetime abstract index and ^, v are the 
tangent-space index. We should notice that A a correspond to the massless U(l) gauge 



field in (|2 . 1|) . The spin connection 1-forms 



[iv J a 



(3.3) 



Here (e^) a form a set of orthogonal normal vector bases. The Dirac equation derived 
from the above action reads 

T a D a ifj - rmp = 0. (3.4) 
Without loss of generality, we choose the orthogonal normal vector bases as follows 

>dr J ' 



'9 



(3.5) 



Then based on the ( |3.3| ), the non- vanishing components of spin connections can be 
obtained as 



Vg TF d r y/g^t(dt) a , 
r g TT d r y/g^(dx l ) a . 



(3.6) 



lr The bulk action of fermion with bulk dipole coupling was discussed in 23 
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In order to remove the spin connection in Dirac equation and to investigate in Fourier 
space, following || , we make a transformation 



~9Q rr )~ 1 e ~ iujt+ik i x ' 



(3.7) 



Considering the rotation symmetry in the spatial directions, we can simply choose the 
momentum along the x-direction k = k x . Then, the Dirac equation can be written as 



^T r d r cf) - i(u + qA t ) V£ ttr V + iksfg^T x (j) - mcf) = (3. 
As in ^ 25 , we will consider the following basis for gamma matrices 

ia 1 ! 



-a 3 l 
-a 3 l 



ia l l 



ii ' 



-a 2 l 
a 2 l 



(3.9) 



where 1 is an identity matrix which have size 2 d 2 3 for d odd and size 2 d 2 4 for d even. 



Now, we set 



then (|3.8|) becomes 



r \ ' ) • 1)10''' ' 



+ qA t )ia 2 <g> ^ J =F ^v 7 ^ 1 ® ( ^ ) ( 3 - 10 ) 



In order to decouple the equation of motion, we set 4>i = ), with / = 1,2. The 
equation of motion yields 



Vi 



{^d r + 2m)6 = [V?*(w + qAt) + (-l) J £V<Fj + [Vs*(w + ?A) - (-l) J AV<Fj £/, 

(3.11) 

where 6" = fj- Substituting the black hole sulution (|2.2|) to the Dirac equation and 
make some algebraic calculation, we have 



(a; + qAt 



-1) J 



+ 



£■ (3-12) 



Asymptotically, when r — > oo, namely towards the boundary, the solution of the Dirac 
equation ( |3.10D can be written as 



'/ -> -V" ( ° ) £//■" 



If we suppose two coefficients in ( |3 . 1 3| ) are related by 



(3.13) 



(3.14) 
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then Green function on the boundary isp4 



G = -iSj°. (3.15) 
So the Green function can be written as 

G= lim r 2m J . (3.16) 



£ 



On the other hand, the ratio G = can be fixed by imposing the in-falling boundary 
condition for <pj at the event horizon, where (pi is behave as 



^ l0C \l) e_1Wr+ ' ( 3 - 17 ) 
with r + = f 4fj. That is to say, this ration can be written as 

G = lim r 2m £/, (3.18) 
which is same as (|3.16|) . The in-falling boundary condition at the horizon yields 

6 |r=r+= i. (3.19) 

4. Numerical calculation 

In this section, we will explore some properties of this system by investigating the 
behavior of the retarded Green function. At first, we try to obtain the fermi momentum 
numerically in the Lifshitz background ( |2.2|) . Before processing, we would like to point 
out that by taking k — > —k, the relation between G\\ and G 2 2 is G^w, k) = Gn(w, —k) 
that can be deduced from the equation of motion ( |3.11Q and the boundary condition 



( |3.19| ). Thus, we only focus on G22 in the following discussion. For convenience of our 
numerical calculation, we can set m = 0, L = 1. 

4.1 Charge dependence 

In condensed-matter physics one more often deals with the situation that under a scale 
transformation the time direction scales twice as fast as the spatial directions, i.e., the 
dynamical exponent is z = 2. So we will discuss z = 2 case with different charge q 
in this subsection. Here, we only study the system with finite temperature. At first, 
we will fix the temperature at T = The behavior of imG22 can be summarized as 
follows: 
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Figure 2: The 3D plot for q = 2, q = 2.5 and q = 3 (w = 1 x 1(T 12 ). 

1. From Figure |1] and Figure 0, we can see that in the limit u — > 0, I mG 22 has a 
sharp quasiparticle peak at the Fermi momentum k F , which agrees with [7|. 
Note that in order to determine k F numerically, we have set u = 1 x 10~ 12 . Figure 
D shows that k F ~ 1.9 for q = 2, k F m 2.5 for g = 2.5 and k F m 3.2 for q = 3. 2 
More accuracy, the sharp quasi-particle-like peak is located at k F = 1.88553916, 
k F = 2.53153814 and k F = 3.18252078 for q = 2, q = 2.5 and q = 3, respectively. 
This implies the Fermi momentum will increase with the increasing of the charge 
q. This result is also consistent with the discussion in 0. Our result also show 
that for smaller charge q < 1, there is no sharp peak at the temperature T = j-. 

2. We present the 3D plot of ImG22{u, k) in Figure 0. It is obvious that the sharp 
quasi-particle-like peak occurs near k F for q = 2, q = 2.5 and q = 3 from the 
three figures, respectively. 

3. Based on the two points above, it is also of interest to explore the behavior of 
ImG22 near the Fermi momentum k± = k — k F . We show the dispersion relation 

u*{k ± ) ~ kl (4.1) 

2 From the third part (q = 3) in Figure 1, we can find there are two peaks in this case. As mentioned 
in l|, when q increased, new branches of Fermi surfaces appear (i.e. when q — 3, the Fermi surface 
has two branches). Another peak is occurred at k ~ —1.8 . But, we only work on the peak located at 
k fa 3.2 in this subsection. 
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Figure 4: The scaling behavior of the height of I mG 22^ *{k±), k±) at the maximum for 
q = 2, q = 2.5 and q = 3. 

in Figure [| In the figure, the blue point is the data of k± and w*. The red line 
is the plot of fitting function of our data. Here we find the exponent a is always 
I by changing q from 2, 2.5 to 3. In other words, The linear dispersion relation 
will be hold for all q in the charged Lifshitz black hole with z = 2. 3 

4. However, to further investigate the dual liquid, we move on to check the relation 
between the maximum height of ImG22 and k± (see Figure |j) 

ImG 22 {cu*(k ± ),k ± ) ~ k~/ (4.2) 

In Figure || the blue point is the data of ImG 22 (u*(kx), k±), and the red line is 
the plot of the fitting function of the data. Interestingly, though kp drastically 
change with q, we find here j3 is around 2 that is almost independent of the charge 
q. This phenomenon is also observed in ||. 

5. In Figure [5], we examine the correctness of a and /3 obtained above by checking 
the scaling behavior near the Fermi surface through fitting the whole curves of 
G%i(l}). For q = 2, q = 2.5 and q = 3, we set k = 1.85, k = 2.52 and k = 3.18, 

3 In the authors found that the dispersion relation exponent a decreases rapidly with increasing 
q. But here, we find that exponent a is f for different values of q. That is to say, the linear dispersion 
relation does not change in the z — 2 charged Lifshitz background when q increases. 



- 8 - 




-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 -0.04 -0,02 0.00 0.02 0.04 -0.02 -0.01 0.00 0.01 0.02 



Figure 5: The plot of ReG22(^>) (blue), ImG2i(w) (orange) and fit function (dashed lines) 
for q = 2 with k = 1.85 (left), q = 2.5 with k = 2.52 (middle) and q = 3 with k = 3.18 (right). 



q 


2 


2.5 


3 


3.5 


4 




1.88553916 


2.53153814 


3.18252078 


3.83819011 


4.49787005 


a 


1 


1 


1 


1 


1 


P 


2.02086 


2.01415 


2.02449 


2.01559 


2.02255 



Table 1: The Fermi momentum and scaling behaviors with different q 



respectively. The fitting function of G22 for both sign of u is 

<fc(«.*J («) 

y ai (-k ± ) 2a ) e 

where a , a\ and 7 are positive constants. The above function has a pole in the 
lower half w-plane at 

u c = - ai (-k ± ) 2a e^ (4.4) 

For Q4.3Q , Re u c gives the location cj*(k±) of the peak and —Ira oj c determines 
the width Y of the peak, 

r = (tan 7 )u;, 2 (£a) (4.5) 

The relation between V and may be the reminiscent of the behavior of Fermi 
Liquid. We emphasize that the function (|4.3|) is the best numerical fits we could 
find and from which we could qualitatively check that the scaling exponent a and 
13 obtained in Figure |] is correct. 

6. Up to now, we have examined the behavior of JmG^ at a low temperature. 
In Figure ^, we compare the spectral function at different temperatures. As 
expected, the high temperature T appears to smooth the peak due to thermal 
fluctuations, which agrees with [|J]. If the temperature is high enough, the sharp 
zero energy peak completely disappears. 

In Table [l], we summarize the results for the Fermi momentum and scaling behav- 
iors. We observe that the Fermi momentum kp increases almost linearly as q increases, 
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limli.,! ■ fj* I 



Figure 6: The plot of JmG 22 H for q = 2(k = 1.87), q = 2.5(k = 2.5) and q = 3(fc = 3.18) 
with T = 4- (solid line), T = (dashed line) and T = Jj- (dotted line). 



while the scaling parameters a and /3 approximately independent of q. According to 
2BI E7[, Landau Fermi liquid is characterized by the scaling behaviors a — — 1. 



However, in our case the exponent (3 ~ 2 although a = 1. From the expression QL5|) , 
it seems that this system characterized by the linear dispersion and the quasi-particle 
width quadratic in frequency resembles the Landau Fermi liquid. This is simply be- 
cause the scaling form is destroyed and Senthil's theory does not apply to our case [^6] . 
Let us explain this point more clearly. It is well known that the standard Fermi Liquid 
form is 

G ^ ^ = 71 tur (4 ' 6) 

UJ — Vp{k — Kf) + Zl 

where the width r ~ u 2 . At first glance, the spectral function yields A(k, u) ~ (k-k F ) 2 ' 
which may give us j3 = 2. However, the proper scaling form for Fermi Liquid is a delta 
function: A(k,u) ~ 5(k — kp — oS) and /3 = 1 for this function. A more important 
parameter to characterize the property of the system is the peak width r U as given in 
fj~5l . Thus, our holographic system preserves two features, i.e., the linear dispersion and 
quadratic quasi-particle width, of Fermi liquid at finite temperature. However, another 
important constraint for Landau Fermi liquid is Luttinger's theorem which states that 
the charge density is equal to the volume enclosed by the Fermi surface. After repeating 
the similar calculation in p8| in which a modified Luttinger theorem is given, we find 



that our system explicitly violates the standard Luttinger's theorem. This is because 
that in the probe limit, the charge density is much larger than the volume enclosed 
by the Fermi surfaces. So our holographic system can not be well described by the 
standard Fermi liquid theory. The similar holographic system of non-Fermi liquid with 
the linear dispersion and the growth of the peak with k± is nonlinear also can be seen in 



25 1 . The Luttinger's theorem can be satisfied when the back- reaction of the fermions 



is taken into account [|17|, |29|, |30 |. 



4.2 Dynamical exponent z dependence 

In this subsection, we study the influence of dynamical exponent z on the behavior of 
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Figure 7: The plot of ImG22(k) for at T = for z = 2(solid line), z = 4(dashed line) and 
z = 6(dotted line). The peaks become smooth as z increases (w = lx 10 2 ). 





the retarded Green function. We know that for fixed r+, the Hawking temperature 
T = j-t z + increases as z increases. As shown in Figure |6|, the high temperature appears 
to smooth the peak due to thermal fluctuations. 

To have a clear picture on the effect of the dynamical exponent z, we would like 
to fix the temperature T instead and let r + vary as z increases. For comparison, we 
will study three cases i.e., z = 2, z = 4 and z — 6. In our numerical calculation, we fix 
q = 3, u = 1 x 1CT 12 and the temperature T = j-. Figure shows that the height of 
ImG22 becomes lower as z is increased. For z — 6, the sharp quasi-particle-like peak 
nearly disappear. 
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Figure 9: The 3D plot of imG 22 0, q) and its density plot for z = 4 at T = i (u = lxHT 12 ). 




Figure |8| shows us that when z = 2 the quasi-particle peaks exist at different value 
of q. Actually, the Fermi momentum kp increases as q increases. The effects of z can be 
seen clearly from Figures [5] and For z = 4 and z = 6 cases, the peak and the Fermi 
momenta are seriously suppressed no matter how big q is. That is to say, in higher 
dimensional spacetime and for larger dynamical exponent z, the exponent z plays a 
dominant role. As z is large enough, the peak will be smoothed out. Meanwhile, we 
also find from the density plot of the three figures that the spacing between the peaks 
increases with z. 
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5. Conclusion 



In this paper, we studied the properties of holographic fermions in charged Lifshitz black 
holes at low temperature for massless fermions. In z = 2 charged Lifshitz background, 
we first solved the Dirac equation numerically for three different cases(i.e. q = 2, q = 2.5 
and q = 3). For q = 2, the sharp quasi-particle-like peak is located at kp = 1.88553916. 
For q = 2.5, Fermi momentum of the sharp Fermi surface kp = 2.53153814. For 
q = 3, Fermi momentum kp = 3.18252078. The value of the kp actually is very 
sensitive to the boundary conditions. When we increase q ,the Fermi momentum kp 
increases with q almost linearly. As q increased, the scaling parameter a remains to be 
1. We clarified that the scaling exponent /3 (/? ~ 2) obtained here does not satisfy the 
scaling theory by Senthil [f26| and the interpretation of /3 should be careful. Further 
numerical computation leads to the conclusion that this system is characterized by the 
linear dispersion and quadratic quasi-particle width which are two important features 
of Landau Fermi liquid theory. However, the system is in fact violating Luttinger's 
theorem. This implies that the holographic system in this charged Lifshitz theory is 
actually fractionalized non- Fermi liquid[|28|]. We then noticed that as the temperature 
increases, the singular zero energy peak is smoothed down due to thermal fluctuations. 
After this, we find the peak of ImG will smooth out when z increases, and the physics 
behind this phenomena need us to further study. It would be interesting to include the 
magnetic field and compare the results of holography with real experiments data in the 
charged Lifshitz backgrounds in the future. 
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